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We provide a general method to find the Hamiltonian ol a linear circuit in the presence of a 
nonlinearity. Focussing on the case of a Josephson junction embedded in a transmission-line res- 
onator, we solve for the normal modes of the system by taking into account exactly the effect of the 
quadratic (i.e. inductive) part of the Josephson potential. The nonlinearity is then found to lead to 
self and cross-Kerr effects, as well as beam-splitter type interactions between modes. By adjusting 
the parameters of the circuit, the Kerr coefficient K can be made to reach values that are weak 
{K < k), strong [K > At) or even very strong [K 2> /t) with respect to the photon-loss rate In 
the latter case, the resonator-|-junction circuit corresponds to an in-line version of the transmon. 
By replacing the single junction by a SQUID, the Kerr coefficient can be tuned in-situ, allowing for 
example the fast generation of Schrodinger cat states of microwave light. Finally, we explore the 
maximal strength of qubit-resonator coupling that can be reached in this setting. 

PACS numbers: 42.50.Pq, 03.67.Lx, 85.25.Cp, 74.78.Na 



> 

m 
(N 
(N 

^' 

O 



X 



I. INTRODUCTION 

With their high quality factors and large zero-point 
electric fields, superconducting transmission-line res- 
onators are versatile tools for the study of quantum me- 
chanical effects in solid-state devices. Resonators have, 
for example, been used to study the strong coupling 
regime of cavity QED in electrical circuits [J [3: to probe 
the displacement of a nanomechanical oscillator close 
to the standard quantum limit Q, to entangle remote 
qubits 0-01 and to implement quantum algorithms d, Q . 

In the same way that Josephson junctions are used in 
qubits for their nonlinearity, these junctions have also 
been used experimentally to make nonlinear resonators. 
With their Kerr-type nonlinearity K{a^ aY /2, these res- 
onators have made possible the realization of Josephson 
bifurcation amplifiers (JBA) to readout the state of su- 
perconducting qubits [lOMlSljOf linear Josephson para- 
metric amplifiers (JPA) |13l - [l7| . of Josephson paramet- 
ric converter (JPC) between microwave photons of differ- 
ent frequencies jl8l - l2C| and the squeezing of microwave 
light [131. these experiments, the Kerr nonlinearity is 
often required to be small with respect to the photon de- 
cay rate k. For example, the nonlinearity is limiting the 
dynamic range of JPAs. Finally, as was theoretically pro- 
posed [2l| and experimentally realized (22j . interrupting 
a resonator with a Josephson junction can also be used 
to reach the ultrastrong coupling regime of circuit QED. 

In this paper, we give a unified description of this sys- 
tem in a wide range of nonlinearity K/ k. Wc first treat 
the very general problem of finding the normal modes of 
a continuous linear circuit in which a Josephson junction 
is embedded. This is done by treating exactly the effect 
of the quadratic (i.e. inductive) part of the Josephson 
potential. The non-linearity is then reintroduced and is 
shown to lead to Kerr nonlinearities and beam-splitter 
type interactions between modes. By adjusting the pa- 



rameters of the circuit, the nonlinearity K can be made 
to reach values that are weak {K < k), strong {K > k) or 
even very strong {K 3> k) with respect to the photon-loss 
rate. In the latter case, the rcsonator-|-junction circuit 
corresponds to an in-line version of the transmon qubit. 
By replacing the junction by a SQUID, the Kerr coeffi- 
cient can be tuned in-situ allowing, for example, the fast 
generation of Schrodinger cat states of microwave light. 
In light of these results, we also revisit the question first 
asked in Ref. j23| concerning the maximal qubit-resonator 
coupling strength that is possible in circuit QED. Wc find 
that the toy model used there and consisting of lumped 
LC circuits current biasing a Josephson junction can only 
be used in a limited range of parameters to predict the 
coupling. 

The paper is organized as follows. In section |lTJ we 
obtain the normal modes of the transmission-line includ- 
ing the effect of the linearized Josephson junction poten- 
tial. We then determine the equivalent lumped-element 
circuit model and reintroduce the Josephson junction's 
nonlinearity. In section IIIIl the three regimes of K/k 
mentioned above are explored. Finally, in section ITVl we 
discuss the question of the maximal qubit-resonator cou- 
pling strength that is possible in circuit QED 



II. HAMILTONIAN OF THE NONLINEAR 
RESONATOR 



In this section, we solve the general problem of a 
Josephson junction embedded in a continuous linear cir- 
cuit. For lumped elements, the approach that we are 
following can be summarized in a few lines. Indeed, as 
illustrated in Fig. [U consider for example a Josephson 
junction of Josephson energy Ej and capacitance Cj in 
parallel with a LC oscillator. The Hamiltonian of this 
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Figure 1. (Color online) An LC circuit in parallel with a 
Josephson junction can be equivalently represented by a LC, 
whose parameters are renormalized by the junction capaci- 
tance Cj and the linear Josephson inductance Lj, in parallel 
with a purely nonlinear element represented by the spider-like 
symbol. 



circuit is simply 



H 



q 

2C" 



2L 



- £;jcos(27rJ/$o), 



(1) 



where C = C + Cj. By Taylor expansion of the cosine 
potential, this can be rewritten in the form 



(2) 



with the renormalized inductance 1/L' = (l/L + 
1/Lj) and where we have introduced 1/Lj — 
{27t/^q)'^Ej the Josephson inductance. In this expres- 
sion, U^S) = £;jS„>,(-l)"+i(2^<5/$o)'"/(2n)! rep- 
resents the Josephson potential excluding the quadratic 
part a d'^. 

As illustrated in Fig. [1] Eq. ^ is simply the Hamilto- 
nian of a linear LC circuit whose parameters L' and C' 
are renormalized by the presence of the junction and con- 
nected in parallel with a purely nonlinear element. Using 
the notation of Ref. 2J], the latter is represented by a 
spider-like symbol. Expressing the conjugate operators 
q and 5 in terms of the oscillator operators a and a\ 
the nonlinear element will contain Kerr- type terms of the 
form (a^a)^ and whose coefficients can be calculated. We 
note that, in principle, the nonlinear terms arising from 
the Josephson potential can be treated to any order as a 
perturbation on the linear Hamiltonian. This approach 
can thus be used in any parameter range. 

In the rest of this section, this procedure will be fol- 
lowed for a transmission-line resonator in which a Joseph- 
son junction is embedded. There, the difficulty comes 
both from the continuous nature of the system and be- 
cause it supports multiple modes. We first give the 
Lagrangian of the resonator including the input/output 
boundary conditions and the presence of the junction. 
The normal modes of the Lagrangian where the Joseph- 
son potential has been linearized are found. Working 
in this basis, the Hamiltonian is then obtained. The 
method presented below is an extension of the work of 
Wallquist et al. [2^ that focused on the first mode of oscil- 
lation of a resonator in the presence of a Josephson junc- 
tion of large Josephson energy (weakly nonlinear regime). 
The reader not interested in the details of the derivation 
can jump to section fll D I which contains the final form of 
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Figure 2. (Color online) a) Discretized representation of a 
transmission-line resonator capacitively coupled to input and 
output ports and with a Josephson junction interrupting the 
center conductor at position xj. b) Lumped element repre- 
sentation of the transmission line -I- junction. The normal 
modes of the resonator, dressed by the junction's capacitance 
and linear inductance, are represented by LC circuits biasing 
the junction's nonlinear inductance (spider-like symbol). 



the quantized Hamiltonian taking into account the non- 
linearity. 



A. Lagrangian formulation and linearization 

Fig. (^k-) presents schematically a transmission-line res- 
onator with a Josephson junction inserted at position xj 
of the center conductor. The resonator, of total length 
2£, is terminated by the input and output capacitors d 
and Co- The Lagrangian of the bare resonator takes the 
standard form 



Cr = 



-i>^{x,t) 



{d^i!{x,t)) 



21 



2L°{x) 



dx, (3) 



where, in the flux representation, ipi^, t) = V{x, t)dt 
where V{x,t) is the voltage, and C^{x) and £° jx) are 
the capacitance and inductance per unit length [2^. For 
simplicity, we will take the parameters of the resonator 
to be piecewise constant. That is, the capacitance C" 
and inductance per unit length to the left (p = I) 
and right (/i = r) of the junction, and the corresponding 
characteristic impedance = i^^/C-^Y^^ group ve- 
locity Vfj, — (L°C|J)~^/^, are not assumed to be identical. 
The generalization to arbitrary C'^{x) and L'^{x) is sim- 
ple 111. 

The input (a — i) and output (a — o) capacitances Cq 
are modelled by 



Co, 



^{Xo,,t)-Vo,{t) 



(4) 
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In this expression, Va (t) is the voltage bias at the port a 
of the resonator and = +i. Finahy, the junc- 

tion is modeled both by its capacitance Cj and Josephson 
energy Ej. In practice, Cj will be a small perturbation 
on the resonator's total capacitance. Dropping this ca- 
pacitance would however be akin to dropping the plasma 
mode of the junction. The contribution of the junction 
to the Lagrangian is then 



Cj 



2 



£;,7cos(27r(5/$o), 



(5) 



where S = ^{x^,t) — ip{xj,t) is the phase bias of the 
junction. The junction can be replaced by a SQUID or 
any other weak link with minimal modification to the 
theory. 

In the same way as for the lumped circuit example 
above, we expand the cosine potential of the junction 
such that the total Lagrangian takes the form 

£ = < + A + £o - (7nl(<5) = /:l - UM6). (6) 

The first term £^ is the resonator Lagrangian includ- 
ing the quadratic contributions of £j and the potential 
Ut^hi^) has been defined as above. We now show how to 
find the normal modes of the linearized Lagrangian £l- 



B. Normal modes decomposition 

In this section, we are interested in finding the orthog- 
onal basis of normal modes of oscillations of the linear 
resonator -|-junction circuit. This is done by solving the 
Euler-Lagrange equation of motions 



dC, 



dC, 



= 0. 



(7) 



In particular, at the resonator ports x = ±£ and the junc- 
tion's position X = xj, Eq. ([7]) determines the boundary 
conditions that strongly influence the resonator mode ba- 
sis. 

Away from the junction and the resonator ports, 
Eq. ([7]) obeys the standard wave equation 



(8) 



whose solutions are left and right movers with the dis- 
persion relation uj^ = fc^w^. Since we are looking for 
modes of the whole resonator, we impose ui = ujr, or 
equivalently that the wavevectors obey Snell-Descartes' 
law krVr = kivi. The field ipi^jt) can be decomposed in 
terms of these traveling modes as 



1p{x,t) = y^^-ipm{t)u,n{x), 



(9) 



onator that we use as a reference and Um{x) the mode 
envelope which we now specify using the boundary con- 
ditions found from the equation of motion Eq. ([7|). 

Indeed, at x = ±£, because of the input and output 
capacitances, the field must satisfy 



4,i-£,t) 



1 



(10) 



da:1p{x,t)\x=e = Vo. 



These equations have homogeneous (Vi^o = 0) and partic- 
ular {Vi^o 7^ 0) solutions. Moreover, the current on either 
side of the junction is equal, which imposes 



(11) 



Without linearization, the last term of this expression 
would contain the full Josephson current /c sin(27r(S/$o)- 

These constraints can be satisfied by choosing 

^ \sm[k^{x + £) - ipl„] ~£<x<Xj 

u,n [X)- ^ j _ ^ ^^^^^ x+<x<£, 

(12) 

with k'j^Vr = kmVi. The normalisation contants Am, rela- 
tive amplitudes Bm, phases f!^ and wavevectors km are 
still to be specified. 

First, answers to the homogeneous solution of 
Eq. (Uni), leading to tan(^f„ = \Za{uJm)/Zf\ with 
Za{uj) = {iujCa)~^- Moreover, the first equality of 
Eq. ([Tl]) fixes the relative amplitude Bm of the modes 
to be 



Bn 



Z° cos[km{xj + £) ~ ^l, 
Zfcos[k'^{xj-£) + ^';^ 



(13) 



This corresponds to the impedance mismatch between 
the two resonator sections at frequency w,„. Finally, an 
eigenvalue equation for the wavevector km is found by 
inserting the envelope function Eq. (|12p in the constraint 
Eq. (ini). This yields 



zl 



tan 



km — {xj - £) + tpm 

Vr 



tan[fc,„(a;j + £) - ipl 



— h f 

(14) 



a transcendantal equation whose solutions are found nu- 
merically. 



with ipmit) oscillating at the mode frequency uJm = kmVi- 
In this expression, km is the wavevector of the left res- 



The remaining parameter to be set is the normalization 
Am- This is done by noting that the Um{x) obey the inner 
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Figure 3. (Color online) Typical example of the first three nor- 
mal mode envelopes of a resonator with a Josephson junction 
at xj = 0. 



proposed in Ref. [2l|, this kink can be used to strongly 
phase bias a qubit and to reach the ultrastrong coupling 
regime of circuit QED. 



C. Hamiltonian of the linearized circuit 

Using the normal mode decomposition Eq. ^ and the 
orthogonality of the mode envelopes Eq. ([T5|) . the La- 
grangian of the linearized circuit can be expressed as 



1 • -0^ • 



This immediately leads to the Hamiltonian 



2Cy 



2Lr. 



(17) 
(18) 



product [27 



dx C°(a;)Wm(x)u„(x) + CiUm{-tjUn{-tj 

1 



(15) 

In this expression, Aum = Um{x~^) — Um,{xj) is an im- 
portant parameter corresponding to the mode amplitude 
difference across the junction. The total capacitance 
Cs = J^^C°{x)dx + d + Co + Cj is here used to 
fix the normalisation constant | A„i | . 

Lastly, from the expressions Eqs. ([T4)) and p5|) . it is 
useful to define the inner product of envelope derivatives 
as they are found to obey a similar orthonormality con- 
dition 



dx 

_,1J(X) 



dxUyn{x)dxUn{x) 



— Au„,Au„ 



(16) 



L„ 



Here, we have defined the mode inductance = Csti;^„ 
corresponding to the effective inductance of the resonator 
mode m taking into account the inductance provided by 
the Josephson junction. 

The first three mode envelopes u„i{x) are shown in 
Fig. [3] for a junction symmetrically located at xj = 0. 
Because they have a finite slope there, the odd modes 
phase bias the junction which results in a kink Aum^odd 7^ 
in the mode envelope. On the other hand, the even 
modes do not feel the presence of the junction at this 
location and are unperturbed, Awm^cvon = 0. As first 



corresponding to a sum of harmonic oscillators. Here, 
Qm = SC/Sipm is the charge conjugate to ipm and Qg^m = 
J2a '^rn{xa)CaVa is a gate chargc associated to the volt- 
age bias at port a. 

As will become clear when reintroducing the nonlin- 
earity, for modes having Aum ^ 0, it is advantageous to 
work with the rescaled conjugate variables (pm = 'ipmAum 
and Pm = Im/Aujn. In this language, the above Hamil- 
tonian takes the form 



[Prn - Pg,r,iY , (f'l 



2CL 



2LL 



(19) 



with = Cs/Au^„ and L'„^ = L„iAu^„ and pg,™ = 
Qg^rn/ Aum- The symbol • restricts the sum to modes that 
are affected by the junction and o to those that are unaf- 
fectedbyit: Ho = E^o(9- - gg.m)V(2Cs) + V'^/(2L„). 
As for the simple lumped-element example presented 
above, the presence of the junction rcnormalizcs the pa- 
rameters and L'^ of the (effective) oscillators. 

To better understand the role played by the junction in 
the different modes m, it is instructive to define the effec- 
tive resonator capacitor Cm and inductor Lm for mode 
m (for which Auj-n 7^ 0) in the following way 



Cm 



' — 



dxC°{x)ul{x) + C,ul{-£) + Coulii), 

9 



'dxUmix)f . 



LO{x 



(20) 



To these two quantities are respectively associated the 
electrostatic and magnetic energy stored only in the res- 
onator and not in the junction. The capacitive and in- 
ductive participation ratio of the junction can then be de- 
fined as 77c,m = Cj/C'm and rji^m = L'm/Lj respectively. 
As it should, from Eqs. ((TS|) and ((T5)) the participation ra- 
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tio of the junction and of the resonator for a given mode 
m sum to unity 



Vc 



Cm 



1. 



(21) 



Both ratios are such that 7yc(i),m — in the hmit where 
the junction becomes a short. Moreover, with Cj <C 
Cs in practice, the participation of the junction to the 
electrostatic energy is small, ?7c,m ~ 0. 

It is important to understand that the participation 
ratio of the junction rji^m can be quite different from one 
mode to the other. As mentioned above in relation to 
Fig. [3J some modes do not feel the presence of the junc- 
tion. Because of this variation in participation ratio, the 
mode frequencies are not uniformly spread, a;,„ 7^ m x wi. 
This inharmonicity can be tuned by changing the position 
of the junction in the resonator or by changing Ej. The 
latter can be done in-situ by replacing the junction by a 
SQUID. As will be discussed below, for some applications, 
this inharmonicity can be advantageous. 



D. Reintroducing tlie nonlinearity 

Now that we have the exact modes of the linearized 
circuits, we reintroduce the nonlinear potential ?7nl(^)- 
In terms of the mode decomposition Eq. (|9]), this takes 
the form 



(-1) 



J>1 



2i\ 



27r 
$0 



2i 



El 



(22) 

Using the rescaled variables defined in Eq. we can 

get rid of the explicit dependance on to write the 

total Hamiltonian in the simple form 



{Pm. - Pg,mY , 4'1 



2L' 



. (23) 

In essence, the system can be modeled by the simple 
lumped-element circuit presented in Fig.[2t). This circuit 
is composed of a discrete set of parallel LC oscillators 
biasing together a purely nonlinear Josephson inductance. 
For the remainder of this Article, we will focus only on 
the modes affected by the junction. 

We now quantize the Hamiltonian and introduce the 
creation (aj^) and annihilation (am) operators of excita- 
tions in mode m: 



0m ~\ 



lhCLuj„ 



{aln + am) , 

1 / t _ 



(24) 



The above Hamiltonian now takes the form 
with 



and 



(-1)'+! /27r 



(25) 
(26) 



i>i 



X El 



2il V*o 



2i 



(27) 



While higher order contributions can easily be taken 
into account, for simplicity, here we assume small phase 
fluctuations and consider only the first contribution to 
the nonlinear term. This yields 



n 



A' 



rnn ■[■ ■[■ 



E 



E 



^Cmmrii^liO^m + (^m*^?^) (28) 
rnn 



where we have neglected terms rotating at frequencies 
faster than jo;™ ~ w„|. The nonlinear terms induce a fre- 
quency shift Lo'^ = u!m ~J2n ^mn On thc modc frequen- 
cies, as well as Kerr cross- Kerr Kmn and beam- 
splitter like C/mn interactions with amplitudes 



Kmm =E'c^mVl,m/fl, 



Clmn =(1 ^ Sim/2) i^K „ijnK mi) 



(29) 
(30) 
(31) 



We express these quantities in terms of the charging en- 
ergy E'(j ^ = e^/(2C^) and the inductive participation 
ratio ?7i^m- It should not come as a surprise that, in the 
same way as for the transmon qubit [2^ , the self-Kerr co- 
efficient (i.e. the anharmonicity) is related to the charging 
energy. Here, however, this nonlinearity is 'diluted' by 
the finite inductance of the transmission line which leads 
to a non-unity participation ratio rji^m of the junction to 
a given mode m. We will come back to the analogies of 
this system to the transmon below. 

Finally, as already pointed out, the Hamiltonian of 
Eq. (pS)) includes only the first contribution cx (5^ of 
the nonlinearity. This expansion is thus valid for weak 
nonlinearities with respect to the mode frequencies, 
Kmm/^m 1- Howcver, as discussed in the next sec- 
tion, this does not prevent the nonlinearity to be strong 
with respect to the photon damping rate k. 
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III. THREE REGIMES OF NONLINEARITY 

In this section, we explore analytically and numeri- 
cally three regimes of nonlinearity. Comparing to the 
resonator photon-loss rate k, we define the weak K < k, 
strong K > K and very strong K ^ k regimes. Here, we 
will choose the junction location along the length of the 
resonator to optimize various quantities. It is interesting 
to point out that this choice is not possible with a A/4- 
typc resonator where the junction is. by default, at the 
end of the resonator. The advantage of this additional 
design flexibility is illustrated in Fig. |31 There, we plot 
the frequency wi and Kerr nonlinearity Kn of the first 
mode of a resonator as a function of junction position xj 
and Josephson energy Ej. The parameters can be found 
in the caption of the Figure. By appropriate choice of 
xj, Ej and, as we will discuss below, the total length 21 
of the resonator, it is possible to reach the three regimes 
mentioned above. 

While the position can only be chosen at fabrication 
time, Ej can be tuned in-situ by replacing the single 
junction with a SQUID with junctions of energy Eji and 
Ej2 and using an external flux ■ In this situation, the 
Josephson potential T-Lj = —Ejcos{2ttS/^o) is replaced 
by lii 



n, = -E. 



cos 



cos 






) sin { 








V *0 y 







dsin 



-£;j($.)cos[27r(,5-,5o)/$o], 



(32) 



where 



Ej{^x) = EjY. cos 




$0 



(33) 



with EjY, = E. 



Ej2, d 



\Eji - Ej2\/Ej^ the asym- 
metry parameter and tan^o ~ rftan(7r$2;/$o)- Below, 
we will drop the phase Sq which can be eliminated by a 
gauge transformation. 

Finally, we note the region of large Kerr nonlinearity 
when the junction is placed close to the of the resonator 
[bottom right corner of Fig. Hb)]. As will be discussed 
in Sect. IIVI in this situation the junction and the small 
segment of resonator to its right essentially behave as 
a transmon qubit, the nonlinearity being related to the 
anharmonicity of the transmon. 



A. K < k: JBAs, JPAs and JPCs 

The regime of weak nonlinearity K < k is realized 
by working at large Ej such that the participation ratio 
rji^rn of the junction, and in turn the Kerr nonlinearity 
Eq. (Unili is small. This regime has been well studied 
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Figure 4. (Color online) a) Frequency lji and b) Kerr nonlin- 
earity Kii of the first mode as a function of xj and Ej. The 
black lines indicate the contour of constant Kh/k = 1 and 
10, with Kxi'ijli^ ~ 1.5 MHz. The dashed white lines indicate 
xj/£ — 0.5, 0.75 and 0.95. These positions will be used below. 
The resonator has a total length 2£ — 1.2 cm, characteris- 
tic impedance = 50 f2 and input and output capacitors 
Ci = Co ~ 10 fF. In the absence of the junction, it has a 
fundamental A/2-mode of oscillation at uj-xj^l'i-''^ = 4.95 GHz. 
Unless otherwise stated, the numerical calculations in this sec- 
tion have been realized with this set of parameters. 



experimentally, starting by the pioneering experiments 
in the late 80's of Yurke et al. with Josephson paramet- 
ric ampliflers (JPA) [H [H. Since, weakly nonlinear 
resonators have been used as a JPAs Il5l-fl7ll . as Joseph- 
son bifurcation amplifler ( JB A) [III Il2l . |24| , as Joseph- 
son parametric converter (JPC) [18l - l20| and to squeeze 
microwave light 
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While in this section we are not 
presenting new ways to exploit this regime, we hope that 
these results, which can be used to predict the important 
high-level system parameters (frequencies, nonlinearities, 
. . . ) from basic device parameters (resonator length, xj, 
. . . ), will prove useful in practice. 

We first briefly comment on the use of this device as 
a JBA, and then move to the JPA and JPC. In the JBA 
mode, the resonator is driven by a tone of frequency 
ujd = A detuned by A from a given mode m. The in- 
harmonicity helps in suppressing population of the other 
modes, for example if |a;„i — uJn=^m\ 7^ 20;^ for the other 
modes n of the resonator. In this situation, the unwanted 
modes can be dropped and the Hamiltonian of Eq. 
takes the simplifled form 



(a„ja,n) 



(34) 



inhere wc have added the drive. Fig. [5] presents the fre- 
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Figure 5. (Color online) a) Frequencies, b) detunings from 
two-photon processes and c) self- Kerr coefficients Kmm of the 
first four normal modes of a nonlinear resonator as a function 
of the Josephson junction position, starting from the center 
{xj = 0) to the right edge {xj = £) of the resonator. Be- 
cause of the shape of the mode envelopes Um{x), moving the 
junction away from the center increases the nonlinearity on 
the higher resonator modes and thus the inharmonicity. The 
dashed vertical line indicates at xj — £/2 corresponds to the 
junction location chosen in Fig. [6] The parameters are given 
in the text. 



qucncics, inharmonicity and Kerr amplitude as a func- 
tion of the junction position (in the range to £) for 
a large junction Ej/h = 636 GHz. As can be seen in 
Fig-EJ)), the detunings to undesirable transitions involv- 
ing two drive photons can be optimized by moving the 
junction along the length of the resonator. 

We now turn to the JPA and JPC modes of operation. 
Here, we will assume that a SQUID rather than a sin- 
gle junction is present. Modulating the flux allows for 
degenerate amplification, as well as non-degenerate am- 
plification and conversion. Indeed, in the presence of a 
small time-dependant flux $rt(t), the Josephson potential 
Eq. ([5^ gains a rf contribution 



EjY. sin((/?^/2) cos(27r(5/$o) 
z [ 

— dcos((^a/2) sin(27r(5/$o) 
where we use the notation — 27r$/$o. Taking iprfit) 



(35) 



(p-[iCos{u)di) this can be rewritten, to second order in 5, 
as 

Hl = -Y. Hg,ne-'"'' + h.c){al + a„0 

m 

- |(5mne-'"''* + h.c)(aj„ -I- a„)(aj, -f- a„), 

ni.n 

(36) 

where the amplitude of the one and two-photon processes 
are 



'27r\ dEjscos{Lpx/2) 



$0/ ^8hC'„,uj„ 
27ry Ejs sm{^,/2) 



= <Prf- 



(37) 
(38) 



For asymmetric Josephson junctions {d ^ 0), both pro- 
cesses are possible. Interestingly, using Eq. ([T6l) . it is 
possible to express gmn in terms of the rate of change of 
the mode frequencies with respect to the external flux as 



(39) 



This relationship, valid for d ^ 1 and $3. 7^ $o/2, was 
experimentally verified in Ref. (20l |. 

We now focus on the symmetric {d ~ 0) case in the 
presence of a non-zero dc component (px 7^ 0. For 
LOd ~ ujj ± Wi, where i, j label two modes of the resonator, 
the Hamiltonian Eq. (|28p including fiux driving can be 
simplified to 



m.,n={i,j} 

h 



(gmne + fi.c)(aj„ + a,„)(4 + 



(40) 



Because of the repeating indices in the sum, the exchange 
rate between modes is gmn- For ~ 2aJm, the last term 
reduces to oc (al^+ai,) corresponding to degenerate para- 
metric amplification |15| . In the non-degenerate mode, 
with the drive frequency « a;^ + a;„ , this term rather 
reduces to cx al^al^ + a„ia„ which can be used for phase 
preserving amplification. Finally, for Ud = a;„i — w„ , we 
find a beam-splitter like interaction aJ^a„H-ama3j between 
modes. Conversion of microwave photons between two 
modes in this JPC mode of operation has already been 
observed [2^. As noted by these authors, the fidelity 
of photon frequency conversion will suffer from the non- 
linearity K„in- Moreover, when operated as a JPA, the 
nonlinearity will limit the dynamic range of the ampli- 
fier and, as discussed in more details in Appendix El the 
number of photons that can be present in the resonator 
before the junction's critical current is reached. The ob- 
jective is thus to increase the coupling gmn while keeping 
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Figure 6. (Color online) a) Frequencies, b) detunings from 
high-order photon processes and c) self-Kerr effects Kram and 
Rabi frequencies gmn for $rf = 0.02$o of the first three nor- 
mal modes of a nonlinear resonator with a large Josephson 
junction, as a function of the external flux The SQUID 
is placed at position xj = ^/2. Parameters are given in the 
text. The vertical dashed line represent at = 0.37$o cor- 
responds to the operating point of the JPC used in Ref. [20l |. 



Kran Small. The added design flexibility provided by the 
choice of the junction's position xj helps in this regard. 
As can be seen in Fig. 21 the dependance of the mode fre- 
quency with respect to Ej (or can be increased by 
moving the junction along the resonator length without 
increasing the Kerr non-linearities significantly. 

With this in mind, we now compare our theoretical 
findings to the experimental setting of Ref. [l^]- There, 
a A/4 resonator, therefore with the SQUID necessarily 
located at one end of the resonator, was used to real- 
ize the JPC Hamiltonian. The necessary inharmonicity 
was realized by varying the characteristic impedance of 
the resonator along its length. In this way, a mode de- 
tuning (cja — W2) — 2(a;2 — wi) = 27r x 240 MHz was 
obtained B By biasing the SQUID at = 0.37$o, the 
JPC coupling was 312 /27r ~ 20 MHz for a flux modula- 
tion amplitude $i.f ~ 0.02$0i while the Kerr coefficients 
were kept relatively small with Kwjl-K ^0.5 MHz and 
is:22/27r ~ 4 MHz. 

Fig. ini presents the same parameters (mode frequen- 
cies w,„, inharmonicity, Kerr coefficients K^m and para- 
metric couplings gran) as a function of the external 



flux In opposition to Ref. ^0|, here we consider 

a A/2 resonator with a symmetric SQUID located at 
xj — 1/2. This location is also indicated by verti- 
cal dashed lines in Figs. S] and [S] We first note that, 
through most of the <^x range, the detunings to unde- 
sired transitions are kept at more than 250 MHz, sim- 
ilar to Ref. [iOl- Moreover, because of the larger flux 
dependance of the mode frequencies, see Eq. for 
the same small rf amplitude — 0.02(f>o, we find JPC 
couphng strengths {gi2, (713, 523}/27r ^ {76,54,86} MHz 
that are about four times as large. Even with these larger 
values, the unwanted Kerr nonlinearity remain small at 
{fs:ii,if22,i^33}/27r ~ {0.21,1.3,0.35} MHz. This in- 
crease in coupling strength over nonlinearity should lead 
to higher fidelities in photon frequency conversion. We 
note that these parameters have only been manually op- 
timized and a more thorough optimization should lead 
to better results. Finally, depending on the flux- noise 
level, it might be more advantageous to work at smaller 
dc flux bias, where the susceptibility to flux noise is re- 
duced, and increase the rf modulation $i.f to keep the 
coupling strength constant. 



B. K > hi: Photon blockade and cat-state 
generation 

In this section, wc will assume that the drive frequency 
is chosen such that the JPA and JPC terms can be 
dropped and that a single mode approximation is valid. 
Similarly to Eq. p4[) . but dropping the mode index m 
and the drive, the Hamiltonian reduces to 



2 ^ ' 



(41) 



^ It is important to note that we take m 
mode, while the authors of Ref. [20II use m 



1 to be the fundamental 
0. 



As already illustrated in Fig. 2] this regime of strong non- 
linearity K > K can easily be reached. Moreover, because 
the Josephson energy is flux-tunable, the Kerr nonlinear- 
ity can itself be tuned. Below, we show how to exploit 
this large nonlinearity, and its tunability, to observe pho- 
ton blockade [s^l and to generate cat states 31 1. 

For the numerical examples, we will use a SQUID of 
total Josephson energy Ejy./2it = 622 GHz interrupting 
the resonator at xj = 3^/4. These parameter are chosen 
such that there is a maximum variation of i^T in a narrow 
range of flux $3,. Indeed, the participation ratio, and 
hence change more rapidly with flux near the end of 
the resonator, as can be seen in Fig. HJd). This narrow 
range should allow for fast-flux tuning to and from op- 
crating points with widely different nonlinearities. How- 
ever, given the dependance of the Kerr effect with the 
Josephson energy, the extremal points are always near 
$2: = and = $0/2, the latter corresponding to 
the point of largest susceptibility to flux noise. To re- 
duce this susceptibility, we will work with an asymmetric 
SQUID of asymmetry parameter d = 5%. In this case, 
the point = $o/2 also corresponds to a sweet-spot for 
flux noise. As can be seen in Fig. [5^), with these param- 
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Figure 7. (Color online) a) Density plot of the average number 
of photons as a function of time and external flux in presence 
of an external drive of amplitude e/27r — 2 MHz and of fre- 
quency LOd ~ uJr — K/2. b) Line cuts for the values of flux 
indicated by I and II in panel a) and corresponding respec- 
tively to K/k ~ 1/25 and ~ 200. 
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Figure 8. (Color online) a) Frequency Ur, radiative relax- 
ation rate k and self-Kerr coefficient K of the first mode as 
a function of external fiux. b) Pulse shape of the external 
flux and Kerr coefficient as a function of time for Schrodinger 
cat state generation, c) Density plot of the numerically ob- 
tained resonator Wigner function immediately after the flux 
pulse acting on an initially coherent state |a = 2). The device 
parameters are given in text. 



eters the Kerr coefRcient can be modulated over four or- 
ders of magnitude, from K{^^ = 0)/27r = 2 x lO'^ MHz 
to K{<^.^ = $o/2)/27r = 20 MHz, over a range of $o/2. 
At its maximum, the participation ratio of the junction 
reaches r]i = 0.6. 

In order to comfortably reach the strong nonlinear 
regime K/k > 1, the coupling capacitances are reduced 
with respect to the previous section to Ci = Co = 2.5 fF. 
This corresponds to a very reasonable photon relaxation 
rate K/27r « 0.1 MHz, or equivalently to = 1/k « 
1.6 fis. With these choices of parameters, the resonator 
can be brought continuously from the linear regime with 
K/k « 1/200 to the strongly nonlinear regime with 
K/k « 200 by increasing the flux threading the SQUID 
loop by half a flux quantum. 

In addition to changing the nonlincarity, flux tuning 
the Josephson energy also changes the mode frequency. 
This is illustrated by the full black in Fig. |5^). Dephas- 
ing due to flux noise due to the finite \duj/d^x\ can be 
evaluated to exceed 10 /us throughout the flux range [28t . 
In practice, the decoherence time should thus be limited 
by relaxation with T2 « 2/k = 3.2 /xs. Finally, the red 
dotted line in Fig. [5^) shows a very weak variation of k 
over the whole flux range. Indeed, while the decay rate 



from the input (left) port decreases with frequency due 
to the Ohmic nature of the bath, the reduction is com- 
pensated by the enhancement of the decay rate from the 
output (right) port as the mode arnplitude at the port 
increases with decreasing frequency [3^ . 

Using these parameters, we show in Fig. [Tji) the re- 
sult of a simulation of the mean photon number (a^a)(t) 
under irradiation with a continuous tone of amplitude 
e/2TT = 2 MHz and frequency uJd{^x) = i^ri^x) — 
K{^x)/2 as a function of time and for various values of 
the external fiux $3,. For small K/k ~ 1/25, correspond- 
ing to line cut labeled I in panels a) and b), the mean pho- 
ton number simply exhibits ringing towards its steady- 
state value. Since the drive is resonant with the |0) O |1) 
Fock state transition frequency, for a large nonlinearity 
K/k ^ 200, corresponding to the line cut labelled II, the 
mean photon number rather shows Rabi oscillations with 
amplitude bounded by (a^a) = 1. To confirm that these 
can be interpreted as Rabi oscillations, the probability 
-Pii^) = is also shown as a blue dashed line in 

Fig. Uh)- This change of behavior from ringing to Rabi 
oscillations is also known as photon blockade, where a 
single photon at a time can enter the resonator because 
of the large photon-photon interaction K (sOj- In this 
regime, the resonator essentially behaves as a qubit with 
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a low anharmonicity K/ujr ~ 0.5%. Leakage to higher 
Fock states can be mmimized by pulse shaping [s^. We 
note that photon blockade was already observed in cir- 
cuit QED using a linear resonator and with a qubit pro- 
viding the nonlinearity, both in the dispersive [sij and 
resonant HH regimes. 

We now turn to cat-state generation. With the res- 
onator initially prepared in a coherent state |a), evolu- 
tion under Hamiltonian Eq. (|4ip for a time r = t: / K will 
generate a superposition of coherent states with opposite 
phases [IH 

l^cat(a)) = - la) + e-"/4|ia)) . (42) 

Here, we suggest to use the tunability of the Kerr coefB- 
cient, with an on/ofF ratio of about 10* for the parameters 
used here, to prepare with high-fidelity this cat state. 

To prepare this state, we first set the external flux to 
~ 0.3<I>o where the resonator is essentially linear. A 
coherent state \a) in the resonator is prepared using a 
tone of frequency uJr for the appropriate amount of time. 
As illustrated in Fig. [5})), a flux pulse '^xit) then mod- 
ulates the nonlinearity K from its small initial value at 
= 0.3$o to its maximum at ^x = 0.5$o, and back. 
By choosing the timing such that the total accumulated 
phase is K{t)dt — tt, the field evolves to a cat state. 
Working with an asymmetric SQUID, the high nonlin- 
earity point where the system spends the most time is a 
fiux sweet-spot, minimizing the effect of flux noise. More- 
over, because of the the large non-linearity K/2'k ^ 20 
MHz, the flux excursion is very short with the required 
phase having been accumulated after a time r = 33 ns 
for the present parameters. And as the resonator is back 
to a small nonlinearity after the protocol, the cat state 
is preserved (up to phase rotations and damping). 

Fig. ISl;) shows the Wigner function of a cat state pre- 
pared using the flux excursion presented in panel b) and 
with an initial coherent state \a = 2). The fldelity 
F[a) = (V'cat(a)|p„um|V'cat(a)) of the numerically ob- 
tained state Pnum to the desired cat state \ipciit{ct)) is 
found to be F(2) « 93.5% for a = 2 and F(^/2) « 97% 
for a = v^j suggesting that photon loss is the main cause 
of error. 



C. K ^ K.: in-line transmon 

By increasing further the participation ratio of the 
junction, the resonator can be made even more nonlin- 
ear. This can easily be achieved by shortening the length 
of the resonator such that its bare fundamental frequency 
(in the absence of the junction) lies comfortably above the 
junction's plasma frequency. In this situation, the first 
mode of the circuit is essentially the junction's plasma 
mode, dressed by the resonator. Focussing on this mode, 
it is useful to write the Hamiltonian in its lumped reprc- 




2i [ml 

Figure 9. (Color online) oji/ojp (full red line) and Ki\/Ec,t 
(green dashed line) vs total resonator length 21. For very 
short lengths, the circuit essentially behaves as a transmon 
of frequency hiOp = ^J9>Ec,tEj and anharmonicity -Ec.t = 
eV[2(2^C''/4-f Cj)]. 



sentation of Eq. (|23p 

n « A.Ecn^ + ^(P^ - E.J cos </>, (43) 

where n — qi/2e is the Cooper-pair number oper- 
ator, conjugate to the phase across the junction. 
The circuit is described by the charging energy Eq = 
(eAMi)^/(2Ci;), the linear inductive energy El = 
(<f>o/27r)^/ (LiAu^), while the Josephson energy Ej is un- 
changed. By taking the one-mode limit and by introduc- 
ing the inductive energy E^ , we have been able to resum 
the cosine potential in Eq. (P^ . In this effective model 
Li/S.u\ plays the role of the inductance of the dressed 
plasma mode. 

With the presence of the inductive 0^ term, Eq. 
correspond to the Hamiltonian of the fluxonium qubit |36l . 
[37| . The analogy should however not be pushed too far: 
the above Hamiltonian is an effective model valid only 
around the frequency of the mode of interest. Indeed, 
while in the fluxonium, a large inductance is shunting 
the junction, there is no such inductive shunt here. As a 
result, the protection against low- frequency charge noise 
provided to the fluxonium by the inductance |37| is not 
present here. In other words, one should be careful about 
reaching conclusions about low-frequency physics with an 
effective model valid only around the (relatively large) 
plasma frequency of the junction. 

As a result, despite the presence of the (j)^ term, 
this system is closer to the transmon [2^ than to the 
fluxonium and, as in Ref. [23| . we will refer to this 
qubit as an in-line transmon. As alluded to before, 
rather than helping, the resonator's inductance dilutes 
the Josephson inductance, reducing the anharmonicity. 
Indeed, the transition frequency for the two lowest ly- 
ing states of the Hamiltonian Eq. can be approxi- 
mated by hujQi w ^J^Ec{El + Ej) and the anharmonic- 
ity a = u!2i " wio ~ Ecrji^i. In comparison, the trans- 
mon transition frequency given by the plasma frequency 
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Figure 10. (Color online) In-line transmon (right) capacitively 
coupled to a low-Q resonator (left). The in-line transmon can 
be coherently controlled, for example, using the right (green) 
control port and the low-Q resonator measured in reflection 
using the left (red) readout port. 



hjjp = \/8EcEj and it's anharmonicity is given by 
a = Ec- For a given charging energy, the anhar- 
monicity of the inline transmon is smaller than that of 
the transmon by the participation ratio rji^i. 

As illustrated in Fig. [3J the participation ratio, and 
thus anharmonicity, is increased by shortening the length 
of the resonator. In the limit where this length goes 
to zero, the first mode has a wavcvector fc — >■ 0, cor- 
responding to an envelope that is constant throughout 
both resonator sections but oscillating out of phase. In 
this situation, the kink at the junction tends towards 
Aui — > 2. The participation ratio rji^i then tends to 
unity, such that El — >■ 0, a — >■ —Ec,t and wi — )- Wp, with 
Ec,T = e^/[2{2eC°/4 + Cj)] the charging energy of an 
equivalent transmon (28j . For this purely plasma mode, 
the in-line transmon reduces to the standard lumped- 
element transmon. 

In the same way as the transmon, the in-line trans- 
mon can be operated in a parameter regime where it 
is protected against charge noise. However, because of 
the finite inductance, the protection here improves with 
the ratio {Ej + El) /Eg rather than Ej j Ec as it does 
for the transmon |37| . In practice however. El <C Ej 
and the extra factor of El should not lead to a signifi- 
cant increase in charge noise protection. However, with 
respect to the transmon, this qubit could benefit from 
lower surface losses. Indeed, with its finger capacitor 
and the associated large electric field, the transmon suf- 
fers from surfaces losses [11]. Here, the capacitive shunt 
is provided by the resonator which does not have (or has 
less) sharp edges and has a smaller surface area. It is 
also possible to further decrease the electric field inten- 
sity by increasing the gap between the center conductor 
of the resonator and the ground planes. With quality 
factors above one million having been achieved with alu- 
minum resonators [39i | and Josephson junctions having 
been demonstrated to be very coherent [11], we can ex- 
pect in-line transmons to have long coherence times. 

Finally, in the same way as the transmon, its in-line 
version can be measured by coupling it to a linear (or 
non-linear [l2|) resonator of typically lower quality fac- 
tor. This is schematically illustrated in Fig. [TUl The read- 
out resonator (left) is capacitively coupled to the short 
in-line transmon (right). The former can be measured 
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Figure 11. (Color online) a) Effective circuit of a junction cou- 
pled to the current fluctuations of an oscillator as suggested by 
Devoret et al. Close to the resonance and with Cj > C 

such that the A/2 mode approximation is valid, this lumped 
circuit is a good representation of the circuit shown in panel 
b). c) Effective circuit for a junction coupled to the voltage 
fluctuations of an oscillator, d) For the physical implemen- 
tation with a transmission-line resonator to be equivalent to 
(c) around resonance, it is required that if £q <^ 2£ for the 
A/2 mode approximation to be valid. As stated in the text, 
the effective circuits of panels (a) and (c) are valid represen- 
tation of the circuits of panels (b) and (d) only around their 
resonance frequency. 

in reflection thought the readout port (left, red) and the 
latter controlled using the control port (right, green). A 
high-Q resonator could be used to mediate entanglement 
between in-line transmons (not shown) (ioj . Finally, we 
note that, using this setup and since the participation 
ratio changes with external flux, the presence of the in- 
ductive term in Eq. could be observed by spectro- 
scopic measurements of the in-line transmon transition 
frequencies with respect flux. 



IV. HOW STRONG CAN THE COUPLING BE? 

The in-line transmon was also suggested in Ref. [11] as 
a way to reach the ultrastrong coupling regime of circuit 
QED, the coupling essentially being between the dressed 
plasma mode of the junction and a dressed mode of the 
resonator. In light of the results obtained in this Article, 
we revisit here this idea. More particularly, we arc inter- 
ested in understanding when the interaction between the 
plasma mode and a resonator mode can be approximated 
by the Rabi-like Hamiltonian 

hK 

Hnahi = hujra''a + hLjpb''b+ — {b^bf + hg{a'' + a){b'^ + b). 

(44) 

where a is an operator of the dressed resonator mode 
and b an operator of the dressed plasma mode. In other 
words, we would like to describe the system using a single 
resonator mode and require the plasma mode to preserve 
a relatively large anharmonicity K. In Ref. [1^, this last 



12 



requirement was not made and we will therefore arrive 
at different conclusions here. 

The circuit analyzed in Ref. (23j is presented in 
Fig. [TTkV As illustrated in panel b), in the absence of 
the junction, this is a lumped element representation of 
a A/2 mode of the resonator with maximums of the volt- 
age at the two ports and a maximum of the current in 
the center. In the presence of the junction, in what limit 
is this lumped-element representation of the continuous 
circuit valid? 

Before answering this question, it is instructive to write 
the Hamiltonian corresponding to this circuit. Using 
the conjugate variables {V'nPr} and illustrated in 

Fig. ITTk). we have 
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(45) 



the three terms corresponding to the resonator, qubit and 
coupling Hamiltonians respectively. In the same way as 
in the previous section, we find that the qubit is renor- 
malized by the resonator's inductance and its Hamilto- 
nian takes the form of Eq. (^5)) with Ec ~ e^/2Cj and 
El = ($o/27r)V4i. 

For Ej,El S> Ec, this in-line transmon is well de- 
scribed by a weakly anharmonic oscillator of plasma 
frequency hiUp ~ y^8Ec{Ej + El) and anharmonicity 
EcEj/{Ej + El)- In this limit, it is useful to in- 
troduce the annihilation (creation) operator fe^''^^ of the 
qubit such that cj) = \/fi/i'2Cjujp){b^ + b). Also writing 

ipr ~ \/ h/ {2CuJr){0'^ + o,), the qubit-resonator coupling 
in Eq. (1451) takes the form 



Hgr = hg{a^ +a){b^ + b). 



(46) 



As in Ref. [Ijl, it is instructive to write the coupling 
strength g in units of the qubit frequency LOp: 
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2ujp V SnaZr 



Ec 



-,1/4 



8{Ej + El) 



(47) 



where Z,. = \/L/C and a ~ Zvac/ {2Rk) is the fine 
structure constant expressed in terms of the vacuum 
impedance Zvac = l/^oc ~ 377 fl and the quantum of 
resistance Rk = h/e^ . Since < ^vaci and given the 
dependence in 1/ ^/a, this scheme appears to allow for a 
coupling ratio g/ujp > 1 comfortably in the ultrastrong 
regime. 

We now discuss the constraints on the circuits for the 
conclusion reached above to be valid. We first note from 
Eq. (j47|) that reaching the ultrastrong coupling regime re- 
quires Ej + El Ec d^l- This is inconsistent with the 



assumption first made when writing the coupling Hamil- 
tonian in the form of Eq. (|46p . More importantly, the 
representation of the resonator in Fig.fTTk). and thus the 
model Hamiltonian Eq. is valid only for a A/2 mode: 
we are here assuming a A/2 mode despite the presence 
of the junction. Indeed, Eq. (|46| assumes only an induc- 
tive coupling between the resonator and the qubit. As 
we have shown above, the resonator mode is dressed by 
the junction which then sees a phase drop proportional 
to Aum- This kink Aum results in a charge coupling 
to the resonator mode while simultaneously reducing the 
inductive coupling Eq. (|T7)) . 

This kink can be minimized if the junction capacitance 
Cj is reasonably large with respect to the resonator ca- 
pacitance 2i!C°, Cj > 2£C°. This minimizes the charge 
coupling to the advantage of the inductive coupling. This 
conclusion is also reached when solving numerically the 
eigenvalue equation Eq. ([T^ requiring ki « 7r/(2^) for 
LUp > uj\/2 [or ^2 ~ 7r/(2^) if ujp < ui\/2]- For exam- 
ple, satisfying Cj > 2^C° with the typical resonators 
parameters given in the beginning of Sect. IIIII requires 
Cj ~ 4 pF. This translates into a very weak charg- 
ing energy Ec/h ^ 5 MHz and correspondingly to a 
small anharmonicity. For the qubit to be in resonance 
with the A/2 mode at ujp/2TT = 5 GHz also requires 
{Ej + El)/Ec ^ 1.25 X 10^. With these numbers, we find 
g/ujp ^ 0.2 from Eq. (j47|) . a value consistent with numer- 
ical calculations. Due to the large junction capacitance 
involved, it can be difficult in practice to reach the ultra- 
strong coupling regime with the setup of Fig. [TTb ) and, 
more generally, realize the Rabi hamiltonian Eq. ([^^ as 
an increase in the coupling is done at the expense of a 
reduction of the anharmonicity. 

To increase the anharmonicity at roughly constant 
g/ojp, a possibility is to slide the qubit away from the 
center of the A/2 mode. In this way, the junction can also 
be phase biased. As illustrated in Fig. [TTli). the largest 
phase bias can be achieved at the end of the A/2 mode 
where the amplitude of the mode envelope is maximal. 
In the limit where the length of the transmission line €g 
on one side of the junction is small Iq ^ 2£, the system 
approximately corresponds to a A/ 2 resonator coupled by 
the junction to an island of capacitance Cs = tqC^ with 
negligible inductance. In this setup, the amplitude of the 
mode envelope at the location of the junction u{t) « y/2 
is only slightly modified by the presence of the junction. 

Starting from the total Lagrangian Eq. dH), neglecting 
the inductance on the right-hand-side of the junction and 
focussing on the A/2 mode, the effective Hamiltonian of 
the circuit is found to be 



H 



Pr 

2CL 



2L' 



27r 

$0 



-^{<t>~^r) 



Cj 



CqC 



jPrq, 
'(48) 

where the effective capacitance and inductance of the res- 
onator are = {21C° + Cj)/2 and l/L'^ = C>r, while 
the qubit capacitance is Cq = Cj+Cg. A lumped element 
representation of this effective Hamiltonian is shown in 
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Fig.fTTbl. In Eq. (|48|) . ipr = V^i^i where lAi is defined in 
Eq. This Hamiltonian essentially corresponds to a 
transmon, of charging energy Ec = e^/(2Cq) and plasma 
frequency huip = ^/SEcEj, coupled to a LC oscillator 
through both phase and charge. In the practical limit 
where Cj ^ Cg <C 2tG^ , the charge interaction is neg- 
ligible. As stated in Sec. IIII C[ the effective circuit of 
Fig.lllb) is a valid representation of Eq. pS)) only around 
the resonance frequency. 

We now evaluate the strength of the dominating cou- 
pling. As done above, we work in the limit Ej/Ec S> 1 
which here allows us to expand the cosine potential of the 
junction. Again, by introducing the creation and annihi- 
lation of the qubit [b) and the resonator (a), we find a 
coupling Hamiltonian of the same form as Eq. (|46p with 

where Z; = yjL'JC, = 2y/Ej(2£C^ is the characteris- 
tic impedance of the renormalized resonator mode. While 
Eq. (|Tf)) was proportional to l/-y/a, here the ratio g/ujp 
is proportional to ^/a. Combined with the small charac- 
teristic impedance Zr < Z^nc, this dependence makes it 
difficult to reach the ultrastrong coupling regime in this 
setting. Indeed, with the same resonator parameters as 
above but now with Ej/Eq ^ 100 and ~ 15 f2, we 
find g/ujp ^ 0.15. However, in contrast to the scheme of 
Fig. II lb ), the anharmonicity here is much higher with 
Ec/h^ 300 MHz. 

For completeness, we now compare the effective model 
of Eq. (|48| with numerical simulations of the full sys- 
tem. Fig. [1^ presents the normal mode frequencies LOm, 
Kerr coefficients Kmm and participation ratios ry;^™ as a 
function of external fiux through a SQUID placed near 
the extremity of a A/2 resonator. The parameters can 
be found in the caption of Fig. [121 In panel a), the full 
lines correspond to the normal mode frequencies ioi and 
UJ2, while the dashed line corresponds to uJr and the dot- 
ted line to uip. These last two quantities are evaluated 
from the effective model Eq. |48|) and, as expected, agree 
well with the numerics away from the region of resonance 
(dashed vertical line). On resonance, an avoided crossing 
is observed and from which we extract g/ujp = 0.12, again 
in good agreement with the effective model. 

Simultaneously, both the non-linear Kerr effect and the 
participation ratio shift from one mode to the other while 
crossing the resonance, denoting the change of charac- 
ter of the excitations from photon to plasmon (and vice 
versa). At = *I'o/2, the participation ratio of the first 
mode reaches its maximum close to unity while the Kerr 
nonlinearity reaches Kn ss Eq /h as expected for a trans- 
mon. The fact that both K and the participation ratio 
are not at their maximum value near = is caused by 
a residual dressing with the second resonator harmonic 
close to 10 GHz (not shown). 

These coupling strengths should be contrasted with 
those obtained using phase-biased fiux qubits [2l| to 




0.2 0.4 0.6 0.8 1 



Figure 12. (Color online) a) Frequencies, b) nonlinear Kerr 
coefficient and c) participation ratios of the first two normal 
modes of the resonator for a junction located near the end of 
the resonator. We have assumed a SQUID of total Josephson 
energy Ejs /27r = 19 GHz and asymmetry d — 5% placed at a 
distance £q = 260 /im from the right extremity of a resonator 
with bare fundamental frequency ijJr/2-K = 4.95 GHz. This 
corresponds to a charging energy of the island of Ec /2n « 
400 MHz [blue dotted line in panel b]. 



lumped [4l| and distributed resonators where ra- 
tios g/ujq ~ 0.12 have been reported experimentally. In 
this setup, it should be possible to comfortably reach 
even larger coupling ratios while maintaining large an- 
harmonicities [2ll. |43|. 



V. CONCLUSION 

We have presented a general approach to find the nor- 
mal modes of a linear circuit in which a Josephson junc- 
tion is embedded. To do so, we included the linear con- 
tribution of the junction as a renormalization of the lin- 
ear circuit parameters. The junction nonlinearity is then 
reintroduced and leads to Kerr-type nonlinearities and 
beam-splitter like interactions between modes. This de- 
scription is most practical for nonlinearities that are weak 
with respect to the mode frequencies, but can still be 
large with respect to the photon damping rate. Indeed, 
we have discussed ways to reach the regimes of weak 
{K < k), strong {K > k) and very strong {K ^ k) 
nonlinearity with respect to damping. These results can 
be used to optimize JBAs, JPAs and JPCs. We have 
also suggested an approach to generate high-fidelity cat 
states by tuning rapidly the nonlinearity of the circuit. 
In the regime of strong nonlinearity, the system behaves 
as an in-line transmon. This qubit could benefit from 
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lower surface losses than the transmon. Finally, we have 
explored the possibility to reach the ultrastrong coupling 
regime of circuit QED with the in-line transmon. 

We became aware of related works while finishing this 
paper [H,!!!]. 
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Appendix A: Critical photon number 

A small Josephson energy typically corresponds to 
large nonlinearities. It however also implies a small crit- 
ical current Ic = 2ttEj/^q, close to which the junc- 
tion will switch to the resistive state. Populating the 
nonlinear resonator with large amplitude fields should 



thus eventually make the junction switch. This max- 
imal amplitude of the field can be phrased, roughly, 
in terms of a critical photon number by requiring that 
{nm.\'HL\nm) < Ej, where |n,n,) is a Fock state of mode 
m. A lower bound for this photon occupation can then 
be expressed as 

As it should, this result is similar to what is found for 
the number of levels ^ y/ Ej / 8Ec in the potential well 
of a transmon qubit. The presence of 77;. ,„ reflecting here 
the reduced participation ratio of the junction. Of course, 
for transmons, exceeding the critical photon number does 
not lead to switching because the junction is voltage and 
not current biased. 

For the parameter regimes studied here, where the ra- 
tio Ej/E'fj ^ ~ 10^ — 10^, the number of accessible states 
is always much larger than unity with n^; ^ 10 in the 
strong nonlinear regime and nj^^ > 10^ in the weakly 
nonlinear regime. 
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